Modeling a nanocantilever based biosensor using a stochastically perturbed harmonic 

oscillator 
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A nanoscale biosensor most suitable for clinical purposes could be based on the use of nanocan- 
tilevers. These microscopic 'diving boards' are coated with binding probes that have an affinity to 
a specific amino acid, enzyme or protein in the body. The binding probes on the silicon cantilever 
attract target particles, and as target biomolecules bind to the cantilever, it changes the vibrating 
frequency of the cantilever. The process of target binding is a random process and hence creates 
fluctuations in the frequency and damping mechanism of the cantilever. The effect of such fluc- 
tuations are given in our model calculations to provide a broader quantitative understanding of 
nanocantilever biosensors. 
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I. INTRODUCTION 

Biosensors are composite devices consisting of biologi- 
cal sensing elements and transducer systems. The work- 
ing principle of these devices involves the binding of the 
desired analyte to the biorecognition element fixed on a 
suitable support matrix connected to a transducer. The 
binding of analytes cause changes in physical or chemical 
properties of bioreceptive elements together with a sup- 
port matrix, which can then be sensed by a transducer to 
generate an electrical signal. This generated signal quan- 
tifies the amount of analyte deposited. Classification of 
biosensors can be based on either their biorecognition 
mechanisms or on the methodology of signal transduc- 
tion4 

Nanobiosensors utilizing nanocantilevers can provide 
extreme sensitivity in the detection of biomolecules (an- 
alytes) down to a single-molecule level. Detecting par- 
ticular biomolecules can help biomedical researchers rec- 
ognize pathogens and diseases during clinical monitoring. 
Like many other detectors, nanoscale biosensors are char- 
acterized by a quantity called the dynamic range. The 
dynamic range is characterized by two extreme points 
starting at the minimum detection limit for detecting 
mass concentration (intrinsic to the detector) to the up- 
per limit when saturation occurs in the detector^ 

Nanocantilever based biosensors are based on the me- 
chanical motion of cantilevers. A cantilever is an oscillat- 
ing mechanical device whose mass continually changes as 
biological analytes attach to it. The attachment of ana- 
lytes leads to a change in the device's resonant frequency. 
The amount of mass deposited on such detectors can be 
estimated by measuring the shift in frequency of the res- 
onator. These detectors possess high sensitivity because 
their intrinsic mass is small but the sticking of analyte 
molecules causes significant mass change and hence pro- 
duces a very large change in frequency. The quality factor 



of such nanodevices is also large and that further adds to 
their sensitivity of detection!^ 

Recently, a slightly different kind of nanocantilever has 
emerged. These nanocantilever biosenesors contain mi- 
crofhhdic channels that are used to detect small mass 
species, e.g., cells, proteins etc. These microfluidic chan- 
nels allow one to operate under a low pressure back- 
ground which enhances the quality factor of detector 4 
In order to measure the change in frequency of nanocan- 
tilevers it is required that system parameters should re- 
main unaltered. Ideally speaking, the analyte once at- 
tached to nanocantilever should not detach or move4 
Such attachment-detachment or adsorption-desorption of 
analyte particles,— their momentum transfer to nanocan- 
tilevers 4 and the random motion over detector 4 causes 
fluctuations in resonant frequency along with shifts in 
frequency. These fluctuations, according to fluctuation- 
dissipation theorem broaden the spectrum of nanocan- 
tilever oscillatory motion. 

In order to develop a theoretical model related to de- 
tection process for these nanocantilevers it is incredibly 
important to keep all macroscopic and microscopic fac- 
tors involved in the detection process in mind. For such 
a tiny system having dimensions on several tens to hun- 
dred nanometers, it is indeed possible to analyze some 
variables macroscopically, but was critical to recognize 
the elements that were less obvious and would be negli- 
gible if one was not working on the nanoscale. Because 
amino acids are one of the main bioreceptors used in 
nanocantilever-based biosensors, the affinity and speci- 
ficity of amino acid sequences was a very important fac- 
tor in our analyses. These amino acids are critically de- 
pendent on ligands to which they are attached. Ligands 
play a major role in the determination of geometric di- 
mensions. This change in geometry could cause receptor 
amino acid sequences to become inactive or even initi- 
ate other events like nonlinear expansion or stochastic 
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motion of the cantilever. These chaotic events need to 
be noted when mathematically modeling nanocantilcver 
biosensors. 

This paper is organized as follows: In section HH we 
give a brief description of how mass deposition causes fre- 
quency shifts in the harmonic oscillator. In section lHTl we 
model nanocantilevers mathematically using a stochas- 
tically perturbed harmonic oscillator. We explain how 
noise introduces measurable random frequency fluctua- 
tions in the oscillating nanocantilever that are propor- 
tional to the mass deposited on its surface. This feature is 
used by biosensors as the detecting mechanism. Finally, 
some concluding remarks are presented in section IIVI 



II. MASS DEPOSITION AND FREQUENCY 
SHIFT 

Nanocantilevers can be simply modeled as damped, 
driven harmonic oscillators. The oscillating mechanism 
we consider in this work consists of a mass-spring sys- 
tem, and its equation of motion will be used to describe 
the classical behavior of nanocantilevers. As it is well 
known, this equation includes a damping and a driving 
terms. Random damping in the oscillator changes its 
frequency. If these frequency changes can be measured, 
then the oscillators (nanocantilevers) can be use to detect 
simple physical phenomena such as mass deposition. For 
a nanocantilever, it is common to measure changes in its 
frequency of resonance. This frequency shift is directly 
proportional to the mass deposited on it as we will see in 
the following. 

For a mass-spring system with mass m and spring con- 
stant k, the characteristic angular frequency under sim- 
ple harmonic motion is given by uiq — ^Jk/m. A small 
change in the mass of the system due to mass deposi- 
tion modifies the frequency of oscillation of the system. 
Taking finite differences of the angular frequency with 
respect to mass and solving for Am, we may get 

Am = Aw - (1) 

For this equation to make physical sense, we assume that 
the elastic and geometric properties of the nanocantilever 
remain unchanged after the small mass deposition. 

When analyte molecules are deposited on a nanocan- 
tilever, there is a momentum exchange between the oscil- 
lator and the particles colliding with it. The process can 
be called a process of random collisions^ These collisions 
bring a random dragging force that affects the nanocan- 
tilever motion. If the concentration of analyte molecules 
is low, then this dragging force can be attributed to a 
molecular drag proportional to the molecular velocity 
and contributing to the dissipative part of the equation 
of motion of the oscillator. This dragging force can be 
written as F ( i = b(dx/db) and corresponds to the damping 
term of the oscillator's equation. When the concentration 



of analyte molecules is low, we can say that Fd is propor- 
tional to the number of analyte particles N a striking the 
nanocantilever per second. However, if the dragging force 
is random, it can be described by the term £(i)(dx/di). 
The term brings random damping to the nanocantilever 
oscillator^ Another possible type of dragging is called 
inertial dragging, and it is proportional to the molecular 
acceleration. We will not discuss this case here. 

Some other type of phenomenon taking place is the 
adsorption-desorption of the analyte molecules on the 
nanocantilever surface^ 7 - Such thermally driven effect 
produces random frequency fluctuations in the oscilla- 
tor^— The process of adsorption-desorption can be mod- 
eled by a molecular flux-dependent adsorption rate and a 
thermally activated rate of desorptionA 7 Another mech- 
anism that produces frequency fluctuations is the one 
occurring with nanocantilevers having microfluidic chan- 
nels. In this case, the fluctuations are produced by the 
diffusion of an adsorbed particle along the microfluidic 
channel inside the vibrating nanocantilever^ All of these 
fluctuations or random stochastic processes mentioned 
above can be incorporated as delta-correlated or expo- 
nentially correlated noise processes (Ornstein-Uhlenbeck 
processes) in the modeling equations. 

III. STOCHASTICALLY PERTURBED 
HARMONIC OSCILLATOR 

Let us consider the simple case of a nanocantilever de- 
scribed by a damped harmonic oscillator of mass m and 
spring constant k, and driven by a sinusoidal external 
force. The differential equation describing the motion is 

d x dx 

m— - + b— — h kx — Fo sm(ut), (2) 
dt z dt 

where Fq and oj are the amplitude and frequency of the 
external driving force, respectively, and b is a positive 
damping coefficient. 

If a system is subject to both random and periodic 
forces, the well-known phenomenon of stochastic reso- 
nance (SR) may emerge.—^ In the following, we will de- 
scribe the effects of including random damping and fre- 
quency as stochastically perturbations in the equation of 
evolution of the nanocantilever. They will enter Eq. ([2]) 
as multiplicative noise such as the new equations of 
motion can be written as 

^+2/3[l + £(t)]^+wfc = Asui(wt) (3) 

for the random damping, and 
x dx 

_+2/3— +^[1 + £(£)]* = A sin^t) (4) 

for random frequency. Here /3 = b/2m is the damping 
parameter, uq — ^Jk/ra is the characteristic angular fre- 
quency in the absence of damping, and A = Fo/m. The 
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random variable will be considered as both Gaussian 
white noise with the correlator 



{£(*)£(*')) - DS(t - 1'), 

and color noise with exponential correlator 



2-\\t-t'\ 



(5) 



(6) 



The parameters D and a represent the white and color 
noise strengths, respectively, and A the correlational de- 
cay rate. For color noise, we are only interested in sym- 
metrical dichotomous noise (random telegraph signal) for 
which £(t) can take one of the values £ = ±a, and the av- 
erage waiting time for each of these states is A -1 . Under 
the following limiting conditions: a 2 — > oo and A — > oo, 
Eq. © reduces to ©, provided a 2 /\ = D (see Ref.[l0|). 



A. Nanocantilever with random damping 

In the absence of an external driving force (A = 0) the 
equation of motion ([3]) takes the form 



~ dx 

D {x} = -2^{t)^ 



(7) 



where the operator On is defined by the following ex- 
pression 



- d 2 nn d 

°^dT 2+ ^dt + ^ 



(8) 



If the random damping is produced by the delta- 
correlated white noise ©, then © reduces to 

^ (x} + 2/3(1 - 2(3D) j t (x) + u 2 (x) = 0. (9) 

Here, we have use the fact that averages split in the 
forroii 
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Equation (j9|) models the collisional damping produced by 
analytes getting stuck on the nanocantilever. If 2/3D < 
1, the presence of white noise will lead to a damping 
decrease (weak noise). On the other hand, if 2/3 -D > 
1, the effective damping turns negative increasing the 
amplitude of oscillation (x) and leading to instabilities 
in the dynamics of the nanocantilever (strong noise) 

When the nanocantilever is driven by an external force, 
a solution to © can be written as 



(x) = B sm(ut + ip). 



(11) 



Upon substituting (fTTj) into ^ , we can solve for the am- 
plitude B and get the following expression 



B = 



A 



[(tu 2 - oj 2 ) 2 + 4p 2 Lo 2 (l - 2(5D) 2 ] 



1/2- 



(12) 



The presence of noise in the system brings fluctua- 
tions in the nanocantilever frequency, resulting in spec- 
tral broadening. Hence it becomes necessary to deter- 
mine the minimum measurable frequency shift that can 
be observed in this noisy environment. The spread in the 
frequency Suq can be obtained by integrating the spectral 
density of frequency fluctuations S(uj) 



Sujn 



S(uj)duj 



1/2 



(13) 



Here, we assume that a measurement of the nanocan- 
tilever frequency was done with a square-shaped transfer 
function over the bandwidth 2Acjo, and centered at ujq. 
Equation (fT3|) is an estimate for any real system. The 
spectral density S (co) is determined by the nature of the 
noise present in the system, and can be determined by 
taking the Fourier transform of the white noise correla- 



tor- 



(£(*)f(0)>e- 



l dt. 



(14) 



At resonance, the second moment of the nanocantilever 
displacement, (x 2 ), satisfies the relation 



(15) 



where k is the Boltzmann constant and T is the absolute 
temperature at resonance. By using (|14p . we can find 
the spectral density corresponding to these displacement 
fluctuations 



S( X 2) - 



const x S(u>) 



[(lu 2 - uj 2 ) 2 + 4/3 2 uj 2 (1 - 2pD)- 



,1/2' 



(16) 



There are two damping mechanisms present in Eq. (|16|) : 
The thermo-mechanical fluctuations governed by the 
damping parameter j3, and the momentum exchanged 
with white noise, whose spectral density is proportional 
to /3-D. In particular, that momentum exchange is re- 
sponsible for taking the system into resonant states. The 
spectral density for the fluctuations in the nanocantilever 
displacements turns out to be a resonance frequency 
curve in the parameter D (the noise strength) which 
shows a maximum at D = (4/3) -1 . This behavior resem- 
bles the stochastic resonance phenomenon (see Fig. [T]). 

In this problem, however, the phenomenon of stochas- 
tic resonance is somehow counterintuitive. The reason 
is that the "unwanted" noise coherently adds up to the 
external force signal, increasing the signal to noise ra- 
tio instead of reducing it.— From Fig. [IJ it is clear that, 
for these nanocantilever based biosensors, the resonance 
peak is not sharp but has a finite width. The peak 
width is determined by the intrinsic damping (/3) of the 
nanocantilever, and the stochastic fluctuations caused by 
analyte molecules sticking on it. To measure the precise 
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shift in the resonance frequency of the nanocantilever, 
the sharpest possible peak from the figure should be the 
appropriate choice. Therefore, knowing precise values of 
(3 and D parameters, and estimating the mass deposition 
correctly becomes critical in designing suitable nanocan- 
tilever based biosensors. 



B. Nanocantilever with random frequency 

In the absence of an external driving force (A = 0) 
equation ((4]) takes the form 

^£ + 2^ +^[1 + ^ = 0. (17) 

Here, the multiplicative noise introduces random 
fluctuations in the nanocantilever frequency. It can be 
shown that, if the noise is white with correlator (|5|), the 
first moment of the oscillator is not affected by the noise 
and (x(t)} = x(t). On the other hand, for the exponen- 
tially correlated noise ([6]), the first moment (a;) satisfies 
the following equation in accordance with the cumulative 
expansion of van Kampeni^ 

^( x ) + (2f3-^ qi )j t (x)+co 2 [l-u Joq2 ](x) = 0, (18) 

where the parameters qi and q^ are given by^ 

qi=2 (£(*)£(* - f'))[l - cos(2uj t')]dt' , (19a) 
Jo 



q 2 = 2 / - t')) sm(2u t')dt' . 



(19b) 



In the limit of white noise both parameters q\ and q 2 
vanish. In other words, white noise does not change the 
width of the resonance profile of the frequency. So the 
rate at which the frequency of the nanocantilever changes 
due to mass deposition remains simple to calculate. 
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FIG. 1. Plot of S/ x 2\ /S(lj) as a function of the driving field 
frequency tj, with parameter conditions: loq = 50 and f3 — 
0.35. Curves A, B, C, D, and E are for D = 5, 10, 20, 30 and 
40, respectively. 



With random frequency, the mathematical expression 
used to find the spectral densities of noise fluctuations 
becomes a bit more complicated than that used in ran- 
dom damping (see Eq. (fT4")) ). Hence, the analytical solu- 
tion (jlip for the first moment in random damping can not 
be used to obtain a solution for random frequency. How- 
ever, for dichotomous color noise, the amplitude associ- 
ated with the first moment shows a stochastic resonance- 
like feature in a and A similar to that observed in (|13p . 
Therefore, a similar analysis can be applied to get the 
resonance profile and find the change in frequency of the 
nanocantilever due to analyte deposition. 

The adsorption-desorption of noise by the nanocan- 
tilever gives rise to fluctuations in its frequency. This 
is caused mostly by the constant bombardment of ana- 
lyte molecules on its surface. This noise mechanism can 
be understood from the following perspective. The an- 
alyte molecules are adsorbed due to their affinity to the 
nanocantilever substrate and are desorbed because of a 
finite temperature change. This creates some fractional 
frequency noise. It is interesting to note that this pro- 
cess of adsorption-desorption of analyte molecules does 
not produce a damping mechanism per se. The ran- 
domness in the sticking and releasing particles on the 
cantilever does not contribute to the average change in 
the energy, but it changes the frequency of the nanocan- 
tilever in a nondeterministic manner.- 7 . Thus, this pro- 
cess introduces a different parametric noise that does not 
culminate into dissipation. 

The nanoscale cantilevers are quite sensitive to the 
adsorption-desorption of noise when compared with con- 
ventional scale cantilevers. The reason is the differ- 
ence between the surface to volume ratio for each type. 
This explains why the number of adsorption locations in 
nanocantilevers is bigger than those on the counterpart. 

The frequency fluctuations caused by noise can also be 
described using some basic considerations different from 
those discussed in (fT5)) . We describe this method briefly 
for the sake of completeness, closely following Refs. 6 
and 0- Let the adsorption rate be R a , which is depen- 
dent on the sticking coefficient of the nanocantilever sur- 
face, and Rd the temperature dependent desorption rate. 
The probability of molecular occupation in a particu- 
lar area is given by p = R a /(R a + Rd), and the cor- 
responding variance in the occupation probability o~ v = 
y/R a Rd/ {Ra+Rd}- The correlation time r c of absorption- 
desorption noise is also given by r c = (R a + Rd)^ 1 ■ The 
spectral density of noise in this case can be written as 



S a (ui) = 



2TTiu^N a a^T c 
[1 + (w - uj ) 2 i 



Am 
m 



— > ( 2 0) 



where Am is the mass of the molecules attached on the 
cantilever surface. The shift in frequency follows equa- 
tion ([T|). The noise variance is a maximum when the 
probability of occupation is 1/2, that is the adsorption 
rate equals the desorption rate. On the other hand, the 
noise is a minimum when the occupation probability is 
either or 1. This noise will be superimposed on the 
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frequency change of the biosensor and hence becoming 
critical in estimating the analyte mass deposition as we 
discuss in the following. 

Integration over the spectral density S a (uj) provides 
some change in the nanocantilever frequency 



Aw 



S a (oj)doj 

-^^[iV a tan-(2.A/r c )] 1 / 2; 



(21) 
(22) 



where A/ defines the width of passband. Hence, the 
change in mass on the nanocantilever is given by 



Ama v [N a tan -1 (2ttA/t c )] 



1/2 



(23) 



Note that this expression for the frequency shift (due 
to mass deposition) may be more realistic than that ob- 
tained in ([l}. 



IV. SUMMARY 

In this work, we presented a realistic model for a 
nanocantilever biosensor based on the description of 
a stochastically perturbed harmonic oscillator. These 
biosensors work by mass sensing the analytes through 
shifts in the characteristic frequency of the oscillator. 
When analytes are deposited on the nanocantilever, they 
bring a variety of noises on this vibrating system which 
gives rise fluctuations in damping as well as frequency 
of the cantilever. The estimation of such fluctuations on 
the spectral response of cantilever is important in finding 
out the exact amount of analyte deposition for clinical di- 
agnostic purposes. In this work we considered frequency 
fluctuations according to different type noise models and 
their effect on the mass deposition on the cantilever. 
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